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We study the problem of closing the detection loophole in three-qubit Bell tests, the experimentally most 
relevant case beyond the usual bipartite scenario, and show that the minimal detection efficiencies required can 
be considerably lowered compared to the two-qubit case. The lowest reported detection efficiency thresholds 
for two and three qubits so far are ~ 66.7% and 60%, respectively. Using the three-qubit W state and a 3- 
setting Bell inequality, we beat these thresholds and with an 8-setting Bell inequality we reach 50.13%. We also 
investigate generic three-qubit states which allow us to attain a detection efficiency of 50% in a 4-setting Bell 
test. We conjecture that the limit of 50% is unbeatable using three-qubit states and any number of measurements. 


I. INTRODUCTION 

One of the most surprising features of quantum mechan¬ 
ics is the prediction that distant parties performing mea¬ 
surements on a shared entangled state are able to gener¬ 
ate correlations which rule out any local hidden variables 
explanation. These nonlocal correlations can be witnessed 
by the violation of Bell inequalities M. By now, many 
Bell experiments using various matter systems have been 
performed (e.g., recently in lUl-d]) providing strong indi¬ 
cation for the existence of nonlocal correlations in nature 
Jit]. However, imperfections in the technical implementa¬ 
tions of these experiments make it possible to reproduce the 
experimental data by local hidden variables model. In or¬ 
der to avoid such a classical explanation, all possible loop¬ 
holes have to be closed simultaneously in a Bell experi¬ 
ment. There are two main technical loopholes, the locality 
loophole and the detection loophole. The former one can 
be closed if there is space-like separation between the ob¬ 
servers such that no signal can propagate from one observer 
to the other. This condition could only be met so far in pho¬ 
tonic experiments JiMl. 

In the present paper, we would like to address the lat¬ 
ter one, the so-called detection loophole. This loophole is 
most relevant in Bell tests which use photons, in which case 
measurements frequently give undetected events. These no¬ 
click events have to be included in the observed data, and 
nonlocal correlations are witnessed detection loophole-free 
only if there is no local hidden variables model of the full 
statistics taking into account the no-click events as well 
iQ. The detection loophole has been closed in differ¬ 
ent physical systems such as ions JH], superconductors Jl, 
atoms a, and more recently in photonic systems as well 

Jill. 

As we have seen, the only system where both primary 
loopholes have been closed are photons, albeit these were 
not closed in the same experiment. Though, important steps 
have been made both experimentally (see references above) 
and both theoretically JT^j, such a loophole-free violation 


of a Bell inequality has not been performed yet. A compre¬ 
hensive review on this subject can be found in Ref. flil . 

Let us mention that closing the detection loophole is also 
relevant from a practical point of view. The more recent 
development of device-independent quantum information 
protocols crucially rely on a detection loophole-free vio¬ 
lation of Bell inequalities. In these protocols, there is no 
need to assume any knowledge regarding the internal work¬ 
ings of the experimental devices used (see iS for a recent 
review of the field). For instance, it would allow two dis¬ 
tant parties to establish a certified secret key in, generate 
genuinely random numbers JTsIl . or perform black-box state 
tomography IB. 

In order to close the detection loophole, we construct 
Bell inequalities which are suited to reveal nonlocality us¬ 
ing detectors with low efficiencies. We will consider the 
relatively unexplored case of three-party Bell inequalities 
involving finite detection efficiencies. In particular, we will 
focus on the case when each party detects particles with 
the same t] detection efficiency. The critical detection ef¬ 
ficiency pcrit, below which nonlocality cannot be guaran¬ 
teed depends both on the Bell inequality considered and the 
quantum state used in the Bell test. 

In the two-party case, rjcrit ~ 66.7% J^ is required to 
violate the Clauser-Horne-Shimony-Holt (CHSH) inequal¬ 
ity JUl] with a partially entangled two-qubit state. For two- 
qubits, to the best of our knowledge, there is no known Bell 
inequality (with possibly more than two settings and more 
than two outputs), which would give a lower threshold. Us¬ 
ing 4-dimensional quantum states and a four-setting Bell 
inequality, this threshold can be slightly lowered (down to 
^ 61.8% 112 ]), however, it is still too high when compared 
to efficiencies achievable with current technology. 

One possible approach to go below these threshold val¬ 
ues is to consider multipartite Bell tests, i.e., more than two 
observers. Buhrman et al. J^ and more recently Ref. J^ 
have showed that an arbitrarily small efficiency 77 can be 
tolerated as the number of parties n and the number of set¬ 
tings m become large. However, these results are interest- 
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ing mainly from a theoretical point of view. Indeed, in the 
experimentally more relevant case of small number of set¬ 
tings, the known results are less promising. For instance, if 
the number of settings per party is fixed to two (m = 2), 
the lowest threshold efficiencies using the Mermin inequal¬ 
ity and its generalized version ll^ were shown to ap¬ 
proach Tjcrit = 50% for large n ll2^ . The same limit can be 
approached if we use the many-site generalization of the 
Clauser-Horne inequality J^. Also, a multipartite two- 
setting Bell test based on single-photon entanglement (i.e., 
a W-state shared between multiple parties) was shown to 
approach rjcrit — 66.7% for large n j^ . These above ex¬ 
amples considered large number of parties and two settings. 
There exist other constructions for class of multipartite two- 
setting i neq ualities (e.g., SHI). Note, however, that due 
to Ref. 11^ the critical efficiency for two-setting inequali¬ 
ties cannot be lower than n/(2n — 1). Hence, none of these 
inequalities may allow us to go below rjcrit = 60% for three 
parties and below rjcrit = 50% for infinite number of par¬ 
ties. 

In contrast to two settings and large number of parties, 
the case of more than two settings per party and moder¬ 
ate number of parties is much less explored. Indeed, in the 
case of three parties (n = 3) and a few number of settings 
TO, which is the experimentally most interesting setup be¬ 
yond the usual two-party scenario, only a few results are 
known. To the best of our knowledge, for three parties the 
lowest detection efficiency is attained in Ref. 11^ giving 
Vcrit = 60% using TO = 2 settings. The aim of this paper 
is to go beyond two measurement settings per party, which 
opens the door to more efficient multisetting Bell inequali¬ 
ties. In particular, we explore numerically the best detection 
efficiencies for the emblematic three-qubit W state iH, and 
also perform detailed numerical search when the underly¬ 
ing state is a more general symmetric 3-qubit pure state. 
Note that the search for critical detection efficiencies using 
the famous Greenberger-Horne-Zeilinger (GHZ) state 
was carried out recently in Ref. attaining the lowest 
efficiency rjcrit = 12/17 ~ 70.59% so far using to = 17 
settings per party (for an explicit construction of the Bell in¬ 
equality, please see the website iHt]). Before this work the 
best bound of rjcrit = 75% for a GHZ state was provided 
by Larsson using the Mermin inequality. 

Here we report a considerable improvement over the 
above values by showing that detection efficiencies as low 
as 50% can be tolerated in tripartite Bell tests featuring a 
reasonable number of measurements. However, our setups 
turn out to be very fragile to noise, hence, we believe that 
the experimental implementation remains a challenging is¬ 
sue. 


II. SETUP 

We consider a Bell scenario with three observers (n = 3), 
Alice, Bob, and Cecil, who carry out experiments in distant 
laboratories. Each observer can choose among to possi¬ 
ble inputs and receive two possible outcomes. Let us iden¬ 
tify the inputs of the three parties with j, fc = 1,..., to 
which correspond to a set of to possible measurements 
{Ai}, {Bj}, {Ck} for each party. Without loss of gen¬ 
erality, we can label with -fl and —1 the two different 
outcomes a, f3, and 7 for the respective parties. The ex¬ 
periment is fully characterized by the conditional proba¬ 
bilities P{af5'y\AiBjCk)- We use the shorthand notation 
P{AiBjCk) = P{lll\AiBjCk) and similarly for a sub¬ 
set of the parties, such as P{AiBj) = P{ll\AiBj) and 
P{Ai) = P{l\Ai), etc. It can be seen that these probabili¬ 
ties fully determine the joint distribution P{aP^\AiBjCk), 
hence it is enough to consider them. 

Throughout this work we stick to symmetric Bell in¬ 
equalities, that is, inequalities which are symmetric for all 
permutations of the parties. In addition, our Bell inequal¬ 
ities will not contain single party marginal terms, they are 
built up only by two-particle and three-particle correlation 
terms. We will also assume without loss of generality that 
the classical bound of the Bell inequalities are zero. The 
Bell inequalities considered in Refs. 1^ . ll2^ are similarly 
restricted. As we will see, this simplification allows us to 
treat the problem with the tools of linear programming. We 
can write such a Bell inequality as: 

m 

^ M^f[P{A,B,) + P{A,C,) + P(S,Q)] 

ij = l 

m 

+ E MglP{AB,Ck)<0, (I) 


where 


M^l=M^=M^l=M^ = Mj^= Mgl 


( 2 ) 


(3) (2) 

and the Bell coefficients Mb/ and Mb are chosen such 
that the classical bound is zero. 


A. Local bound 

Let us first compute the local limit of the above Bell in¬ 
equality ([T]i allowing any classical mechanism. In order 
to do that, it is enough to consider deterministic strategies: 
Each of the parameters at, hi and Ci, where i runs from I to 
TO, may take the value of either 0 or 1, and a deterministic 
strategy is defined by a particular choice. This corresponds 
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to a definite outcome for each measurement value for each 
party. For example, = 1 means that the probability for 
Alice to get the value +1 for her fth measurement is one, 
that is P{Ai) = 1. 

To set the classical bound of the Bell inequality ([T]) to 
zero, we must ensure that 

m m 

'y ' + aiCj + biCj) + 'y ^ M-f.aibjCk < 0 

i,j—l i,j,k—l 

(3) 

for all deterministic strategies. Eq. Q gives 2^"* linear con¬ 
straints for the Bell coefficients. Due to the permutational 
symmetry in Eq. (|2]), two strategies which may be derived 
from each other by swapping the strategies of any two par¬ 
ticipants (e.g. by swapping the values of and bi) lead 
to the same constraint, which makes it possible to reduce 
the number of constraints. Also, Eq. (O is trivially fulfilled 
for any strategy assigning nonzero values for only one of 
the participants. We note that it follows from Eq. Q that 
Mk ^ < 0. We may get this from strategy = l,bj = 1, 
while all other a and 6, and all c values are zero. 


B. Quantum bound 

Now let us consider the quantum case. The maximum 
quantum violation of a two-outcome Bell inequality (i.e. 
the one presented in Q) is always attained by von Neu¬ 
mann measurements IITOll . Moreover, it is sufficient to re¬ 
strict ourselves to pure states that is p = |'0)('!/;|. Then 

P{A^Bj) = {tp\Ai 0 Bj 0 J|V') 

P{ACj) = {2jj\Ai ® i ® Cjl^p) 

P[BiCj) = ® Bi® Cj'IV') 

P{A,B,Ck) = {ii^\A,®Bj®Ck\'<p), (4) 

where Ai, Bj and Ck are the measurement operators of Al¬ 
ice, Bob and Cecil, respectively, projecting onto the sub¬ 
space corresponding to outcome -|-1 in the subspace of the 
participant concerned, and I is the unity operator in that 
subspace. Along this study, we will restrict ourselves to 3- 
qubit states, hence the measurement operators Ai, Bj and 
Ck are in fact projectors in the qubit space. 

C. Quantum case with limited detection efficiency 

Let us consider the quantum case when all participants 
detect their particles with the same limited detection effi¬ 
ciency rj. As a side remark, we note that interesting results 
have been obtained in the asymmetric case, that is, when the 


parties feature different efficiencies EH or when measure¬ 
ments corresponding to the same party have different effi¬ 
ciencies EH- In our symmetric scenario, the participants 
agree to output —1 in case of no detection. In this case, we 
get the joint probabilities of detecting outcome -1-1 by two 
and by all three participants if we multiply the probabilities 
of Eq. (IHi by and by rj^, respectively, that is by the prob¬ 
ability of the detection of the particles concerned. Then the 
condition for the violation of the Bell inequality in Eq. O 
can be written as: 

(V'lTWJV') = > 0, (5) 

where 

m 

^ M^f 0 Bj 0 /|^/:) 

-I- {'4>\Ai 0/0 C'jlV') + (V'l-i® Bj 0 CfclV')) 

( 6 ) 

m 

^ M\%{^P\A®Bj®Ck\i^), (7) 

and AArj is the effective Bell operator at p efficiency. As 
we have shown earlier, < 0, therefore < 0. 

Therefore, if rj is very small, according to Eq. ([SJl, there 
is no Bell violation. The critical detector efficiency, above 
which the violation may be detected is: 

To find the Bell inequality which minimizes ijcrit in case 
of a particular choice of the state and the measurement op¬ 
erators is a problem of standard linear programming. To 
ensure that the classical bound is zero, the set of linear con¬ 
straints given by Eq. ([T]) must be satisfied. As the Bell co¬ 
efficients may be multiplied by any positive number, we 
may fix the norm by fixing the value of M.'Ai. We may 
choose any negative number. In particular, let us choose 
= —1. This provides an additional linear constraint. 
Then we must maximize which is a linear expres¬ 

sion for the Bell coefficients. The symmetries according to 
Eq. (|2]i are further linear constraints to be enforced, but in¬ 
stead of doing that, we may restrict ourselves to coefficients 

(3) (2) 

i ^ 3 Cl k and Mlj , with i < j, and rewrite 

the constraints and the expression to be maximized in terms 

of these independent parameters. This way we get a much 

smaller problem to solve. 

Let the set of measurement operators be the same for all 
parties, and let us confine ourselves to real measurement 
operators. This particular restriction was also proved to be 
useful in other studies for exploring nonlocality of the W 
state EH . In this case the operator Ai = Bi = Ci can be 
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characterized by a single real variable 

^i|0) = ^(1 - cos4>i)|0) - isin$i|l) = c“|0) + Si|l) 

^i|l) = sin$i|0) + i(l + cos4>i)|l) = Si|0) + c+|l). 

(9) 

If = 0, the measurement gives value +1 with probability 
one for the |1) state. 

Let the quantum state be also symmetric in terms of the 
permutations of the parties. One such a state is the 3-qubit 
GHZ state U^l, which case has been already investigated 
thoroughly HH]. In this paper our primary concern is 
the 3-qubit W state ll^ but we also study generic symmet¬ 
ric 3-qubit states. In the following section we focus on the 
W state (Sec. Hilt and then we move on to investigate the 
more general case in Sec.|IV] Our main results concerning 
the found detection efficiency thresholds are summarized in 
Table I and Table II for the W state and the generic 3-qubit 
states, respectively. 

III. DETECTION EEEICIENCIES USING THE W STATE 
The W state is defined by ll^ : 

|IL) = ^(1001) + 1010)-f |100)), (10) 

where we have used the shorthand notation: 

|q;/ 37 ) = |a) (g) 1/3) (g) I 7 ). ( 11 ) 

Now, by using Eqs. (19110111b . it is straightforward to cal¬ 
culate the quantum conditional probabilities appearing in 
Eqs. (I6l7]i : 


will certainly get the Bell inequality that gives the smallest 
critical efficiency with the \W) state. Then for the known 
inequality we may calculate the optimum measurement an¬ 
gles. Also, due to the tightness, the Bell coefficients can 
always be normalized such that they are integer numbers. 

Next we list our results for different number of settings, 
where the numerical study was carried out up to 8 settings 
per party. 


A. W state, m = 2 

Eor two measurement settings per party we got the fol¬ 
lowing Bell coefficients: 

= -1 = 2 m[% = 1 = -1. (13) 

Here we only show the values of the independent Bell coef¬ 
ficients, that is with i < j and with i < j < k. 

The values of the coefficients that can not be derived from 
the coefficients given above by some permutation of the par- 

f 2) 

ties (e. g. M {2 ) are zero. This inequality is equivalent to the 
inequality 22 in the list of Sliwa 14311 . The optimum angles 
for this inequality are $1 = 2.28059 and <I >2 = 0.33432, 
and the critical efficiency is rjcrit = 0.83747. We will see 
later that the \W) state is not the best choice for this in¬ 
equality. 


B. W state, m = 3, 4, 5 

Eor TO = 3 and to = 4 we have got the inequalities with 
the smallest rjcrit if we have chosen the measurement angles 
small. In the case of to = 3, the nonzero independent Bell 
coefficients of this inequality are: 


{W\Ai g Aj g i\W) = i(2siSj + G c+ -f c+c^- -f c^ Cj ) 
2 

{W\Ai g Aj g Ak\W) = - {c~SjSk + s^cjsk + SiSjC^) 

+ \ + C“c+c^ -I- cfcjcl). (12) 

If the number of measurement settings per party is small, we 
can scan the space of measurement angles with an even step 
size, and solve the linear programming problem for each set 
of angles. In each case the optimal Bell inequality we ar¬ 
rive at has to be a tight one in the symmetrized probability 
space. We refer to Ref. ll^ for the framework of sym¬ 
metric Bell inequalities and to further studies which makes 
use of this framework reducing considerably the com¬ 
plexity of the problem. There is a finite number of such 
inequalities, so we get the same solution for a whole range 
of angles. Therefore, if our step size is not too large, we 


= -6 M^f = -3 
-^233 = 2 , 


-^123 “ 3 


M 


(3) 


223 


= 2 

(14) 


while for to = 4 we have got: 



= -6 


= -2 

iV(Zli2 - 

= 6 

-'“114 — 


= 6 


= 3 

IVl 124 - 

= 3 

■'“134 — 


= -6 

CO CO 

= -1 

■^^^334 " 

= 2 

■'“344 ~ 


Eor both inequalities the optimal angles for all measurement 
settings approach zero near the threshold efficiency. This 
observation allows us to make some analytical considera¬ 
tions. 

Let X be small, and let us consider the measurement an¬ 
gles proportional to this small number, that is = (j)iX. 
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Then, Eq. (fT2l i may be approximated as: 


{W\Ai ® 0 i\w) ^ -[1 - cosa;((/)i + 

(16) 


0 Aj 0 AfclhE) ss 4>i‘Pk + 4^j4^k)'^■ 


(17) 


We have neglected terms sixth and higher order in x. 
We have used Eq. @ defining the quantities appearing in 
Eq. (IT 2 I 1 . which may be approximated at leading order as 
Si Ri —(pixj2, fv 1 and c~ « (fPx^/A. Also, it is easy to 
see that 2siSj+c“c^ + clc~ = [1 — cos(<f>i + $j)]/2. Due 
to Eq. (fTTI) . Ad (see Eq. (I7]i) is fourth order in x. Then, 
according to Eq. ([Sll, we may only get a finite value for rjcrit, 
if Ad(^^ defined in Eq. (| 6 ]l is also fourth order in x. This is 
true if whenever the Mb ' Bell coefficient is not zero, the 
corresponding measurement angles satisfy (j)i + 4>j = 0 . 

We may get the Bell inequalities of Eqs. (I14I151 I by solv¬ 
ing the linear programming problem using the small an¬ 
gles limit, that is Eqs. (116117b . when calculating and 
and dropping the overall factor x^. In case of to = 3 
(Eq. fT4l) . we take (j>i = 0 and (j )2 = —(t >3 = 1 (that is we 
choose X = $ 2 ). This way, there are no free parameters 
( 2 ) ( 2 ) 

left. With this choice M^\ and M 23 may take a nonzero 
value, as$i = —$1 = 0 and $2 = — $3 = tc. Indeed, 
these are the nonzero coefficients in Eq. (fT4b . Actu¬ 
ally, the solution of the linear programming problem in this 
case is not unique, there are other Bell inequalities leading 
to the same rjcrit- We have shown the one having the small¬ 
est number of nonzero Bell coefficients. Now, from Eq. ([^ 
we can easily calculate the value of rjcrit- Using the mea¬ 
surement angles defined above, (lU|Ai 0 Ai 0 I\W) ~ 0 
and (W\A 2 0 A 3 0 i\W) ~ a:^/48 (see Eq. [T^. Eur- 
thermore, (lU|Ai 0 A 2 0 A 3 IIU) « {W\A 2 0 A 2 0 
iglVE) = {W\A 2 0 i 3 ® h\W) ^ x^/A&. Also, due 
to the permutational symmetry of the state \W), the ma¬ 
trix elements are the same for all permutations of the op¬ 
erators. Therefore, by substituting the values for the mea¬ 
surement angles and the Bell coefficients into Eq. (| 6 ]l, we 
get = — 18x^/48. Similarly, from Eq. (|7]i, we arrive 
at AdU) = 30a;^/48. Therefore, rjcrit = 3/5 = 0.6. We 
have noted that this is not the only Bell inequality with the 
same threshold efficiency. The reason is that the quantum 
value does not depend on M^^/ *^60 matrix ele¬ 

ments they are multiplied with are zero being $1 = 0. The 
requirement of zero classical value does not define uniquely 
these coefficients. 

In the case of to = 4, similarly to Eq. (El) for TO = 3, 
Eq. (ITSl) can also be derived by using the small angles limit. 
Now, we choose the measurement angles $1 = —^2 = x 
and 4)3 = — 4)4 = Xx. Now we have a single parameter 
A. It is enough to consider |A| < 1. We get the required 


inequality if we choose any value for A between 0.21 and 
0.78. We show in the Appendix that the optimum is A = 
0.466715, which is a root of a fifth order equation, and then 
Vcrit = 0.509036. 

We have also derived the optimal m = 5 Bell inequality 
similarly to the smaller ones in section HIIB I with measure¬ 
ment angles 4)i = 0, 4)2 = — 4)3 = x and 4)4 = — 4)5 = 
Xx. It turned out to be equivalent to the to = 4 case, so we 
got no improvement on the critical efficiency. 


C. W state, rriA = 3 and ms = me = 2 


If we do not require permutational symmetry, we may 
create a Bell inequality with the same rjcrit = 0-6 as for the 
rriA = rriB = rnc = 3 case using the \ W) state with only 
two measurement settings for Bob and Cecil. As before, 
Alice’s measurement settings Ai, A 2 and A 3 are character¬ 
ized by 4)/^ = 0, ^2 = X and 4)^ = —x, respectively. 
However, for Bob and Cecil the measurement angles will 
be chosen as 4)f = 0, 4)f = x and 4)f = 0, 4)^ = —x, 
respectively. Here we used the upper indices to distinguish 
between the parties. The asymmetric inequality will have 
the same quantum value as the symmetric one for any rj, 
if the sum of the Bell coefficients multiplying matrix ele¬ 
ments that have the same numerical value are the same for 
both inequalities. At the same time we must ensure that 
the classical bound is also the same, that is zero. In the 
case of a known symmetric inequality, these requirements 
define a set of linear constraints for the coefficients of the 
asymmetric one. It is a problem of linear programming to 
decide whether these constraints can be satisfied or not. In 
the present case the problem is solvable, the simplest Bell 
inequality we have got, after dividing each coefficient by a 
factor of six is: 

S = -P(ll|AiHi) - P(ll|AiCi) - P(ll|PiCi) 
-P(11|A3P2) - P(11|A2C2) - P(11|P2U2) 
-hP(lll|AiP 2 C' 2 ) -h P( 111 |A 2 Pi(A 2 ) + P(111|A3P2C'i) 
-fP(lll|A2P2C'2) + P(111|A3P2C'2) < 0. (18) 

We have also tried to derive asymmetric Bell inequalities 
with a smaller number of measurement settings for some of 
the parties from the to = 4 case given by Eq. (Ell, and also 
from the inequalities we will show later, but we have found 
no solution for the problem involved. 


D. W state, m > 6 

Eor TO = 6, using 4 )i = —^2 = x, ^3 = —4)4 = 
jix and 4)5 = —4)6 = irx, we got a new inequality, with 
Verit = 0.502417, marginally better than before. Now the 
optimal choice for the parameters is p = 0.495815 and iz = 
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0.295435 (see Appendix). The nonzero independent Bell 
coefficients of this m = 6 inequality are: 



= -18 

JW 34 

= -18 

-''^56 

= -18 

Mill = 18 


= -18 

-'’^122 

= 18 


= 9 

Mill = 9 

-^'^-^136 

= -9 

-''^^156 

= 8 


= -18 

Mill = -9 

-^^-^256 

= 4 

-'^^334 

= 18 

-'^^336 

= -18 

= 18 

■^^'^345 

= 9 

-'^^346 

= 9 

-'^^356 

= 1 

Mi% = -18 

■^^^456 

= 5 

-'^^556 

= 4 

-'^^566 

= 8 , 

(19) 


For TO = 7 we have got no further improvement. 

For TO = 8 there are three free parameters. The an¬ 
gles are given as = —^2 = x, ^3 = —$4 = px, 
4)5 = — 4>6 = ax and $7 = —<I)g = tx . The optimal 
choice of the parameters is p = 0.498442, a = 0.306395 
and T = 0.169989. Then pcrit = 0.501338. The nonzero 
independent coefficients are: 


Mil'> = -6 

1 

II 

-'"56 

= -6 

II 

c7 00 

Mill = 6 

Mill = -6 

Mill 

= 6 

Mill = 

Mill = 3 

Misl = — 3 

Mlf^ 

= -6 

mI% = 

Mill = 6 

M^^^ — —6 

-'"336 “ ° 

Mill 

= 6 

-'"345 ~ 

mI% = 3 

— —3 

-'"358 “ 

Mill 

= 2 

-'"445 ~ 

-^‘^467 ~ 

— 2 
-'"478 “ ^ 

-'"556 

= 6 

-'"558 — 

— 6 
-^^^566 ~ ^ 

Ar(3) _ 0 

^“567 — ''5 

-'"568 

= 3 

M^l = 

A/r(3) _ 

— 1 

-'"678 — ^ 

-'"778 

= 2 

Mill = 


( 20 ) 

We have not tried any larger numbers of settings, the 
number of constraints are too large. We may have got fur¬ 
ther improvement, but we do not expect we could go below 
0.5 with the critical efficiency. 

We summarized critical detection efficiencies we found 
in this paper for the 3-qubit W state in Table|I] 


IV. DETECTION EEEICIENCIES EOR SYMMETRIC 
3-QUBIT STATES 

By considering a more general symmetric state we have 
been able to reach pcrit = 0.5 exactly already with to = 4. 
But we found improvement even for to = 3. The state 
considered is: 

I')/') = cosa|VF)-I-sinajlll). (21) 

This state is also symmetric for the permutations of the par¬ 
ties, therefore, the matrix elements of the tensor products of 


settings 


equation 

222 

0.83747 

list 

223 

0.6 

GU 

333 

0.6 

d 

444 

0.509036 

d 

666 

0.502417 

d 

888 

0.501338 

d 


TABLE I. Table for critical detection efficiencies using the W 
state. The numbers in brackets refer to the Bell inequalities. In 
the first column ijk refers to the respective number of settings i, 
j, and k for the parties Alice, Bob, and Cecil. Detection efficiency 
thresholds are rounded up to five digits. 


single party operators will not depend on the order of those 
operators. In the above state we find that the weight of | 111) 
goes to zero as the threshold efficiency is approached. Like 
before, the measurement angles also vanish at porit- 

Now, besides the matrix elements calculated with the 
I IF) state (see Eq. (fT2l i. and Eqs. (I16ll71 i) for the condi¬ 
tional probabilities of Eq. (|4]i appearing in Eqs. (I6l7l i we 
also need: 

- - - 1 

{W\Ai 0 Aj 0 /| 111 ) = ~ 4 ^'^*'^'^ 

x^ 

{W\Ai 0 Aj 0 Afe|lll) = -^^{cfsjSk + SiC+Sfc -I- s^SJcl) 

~ 4>j4>k) 

(lll|i* 0 Aj 0 /|111) = C+C+ « 1 

(lll|i, 0 A, 0 ifcllll) = C+C+C+ 1. (22) 

These matrix elements, including their limits for small 
angles, may be calculated similarly to the ones given in 
Eqs. (II 2 II 6 II 7 I 1 . We have also used the same notations. 
Eor small angles we have kept only the leading order terms. 

Eq. (I 22 I) shows that the matrix elements (IF|A4,;|111) and 
(lll|Ad^|lll) of the effective Bell operator (see Eq. (|5]l) 
for small measurement angles, that is for small x, will be 
second and zeroth order in x, respectively. With |'(/;) given 
in Eq. (|2T]) we can write: 

= COS^ a{W\Mr,\W) 

+ 2sinacosa(IF|A4,;|lll) + sin^ a(lll|Ad^|lll). 

(23) 

Let us make the same restriction as before, namely let 
( 2 ) 

Mlj = 0, whenever the corresponding measurement an¬ 
gles do not satisfy ({>1 + ({>2 = 1. which makes sure that 
(w\JAr^\W) fourth order in x. Then ('0|Ad,,|'0) is also 
fourth order, if the mixing angle a is taken proportional 
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with x^, that is a = ax^. For small x we may write; 


ft! (VF|7W^|VF) + 2x^{W\M^\lll)a 
+ x‘^{lll\Mr,\lll)a^. (24) 


With this choice, all matrix elements appearing in and 
Al(3) according to Eqs. (1617b are fourth order in x, and we 
may derive the Bell inequalities with the smallest critical 
efficiency using linear programming exactly the same way 
as we have done with the \ W) state. There is one extra pa¬ 
rameter a characterizing the mixing angle. From Eq. (l24l) it 
is easy to determine the optimum choice for this parameter. 
The equation defines a parabola as a function of a, and its 
maximum value is given as 

a = -(VF|7W^|lll)/a;2(lll|7W^|lll). (25) 

We note that < 0 for small a;, which follows 

from the condition that the classical bound is zero, and that 
all values of matrix elements involved are approximately 
one (see Eq. (l22li '). Then the quantum value with the opti¬ 
mum a may be written as; 


ii’lMnl'tp) 


{W\Mr,\W} 


{W\Mr,\ni) 


(26) 


The optimum value of a depends on the Bell coefficients 
to be determined, so what we can do is to try some initial 
values for a, determine the Bell inequality with linear pro¬ 
gramming, calculate the optimum a for this inequality, then 
repeat these steps until convergency, which typically means 
just a few iterations. 

Let us first start with the smallest number of settings con¬ 
sidered; 


A. Symmetric state, m = 2 

Choosing the parameter a according to (|25]) in the state 
Itjj) in Eq. (ITTI i. noting that a = ax^, we may get rjcrit = 
0.6 in the limit of small measurement angles with m = 2 
measurement settings per party. If we choose 4>i = 0 and 
$2 = X, we get the same Bell inequality as we got with 
the \ W) state, we have shown in Eq. (fTJb . The marginally 
small admixture of the |111) state lowered the value of rjcrit 
from 0.83747 to 0.6, with considerably different measure¬ 
ment angles. The inequality is the same as the three party 
one given by Larsson et al. and which is number 22 
on the list of Sliwa j^ . However, in lEsIl the state they 
considered is the |000) state with a very small admixture 
of the \ W) state, that is their state approaches a separable 
state at the threshold efficiency. Also, in their case, the sec¬ 
ond measurement angle is zero, and not the first one. Sur¬ 
prisingly, their very different solution does lead to the same 



FIG. 1. The maximum violation as a function of the detector effi¬ 
ciency for the Bell inequality with two settings per party. 


Vcrit = 0.6. We have calculated the maximum violation 
of the inequality numerically for several detector efficien¬ 
cies above rjcrit- It turned out that it is always enough to 
consider permutationally symmetric real states and to take 
the same real measurement operators for each party. There¬ 
fore, the state can be written as a linear combination of 
\W), I 111) and |000) (the fourth independent real symmet¬ 
ric state can always be eliminated by an appropriate choice 
of the local coordinates). 

The maximum violation as a function of the detector 
efficiency is shown in Eig. (ID- Near the threshold effi¬ 
ciency the maximum violation scales as the third power of 
Ap = rj — rjcrit- The optimum state approaches the \W) 
state, while the coefficients of the |111) and the |000) states 
are proportional to Ary and Ary^'^^, respectively. If we take 
the coefficient of the |000) state exactly zero, the maximum 
violation remains basically the same. Near the threshold the 
difference is negligible, and it is just a little more than 3% 
around rj = 0.9. Therefore, the optimum solution may be 
reproduced almost exactly with the state we have consid¬ 
ered in the present paper. Near rjcrit the measurement an¬ 
gles T*! and $2 scale as Arj^^'^ and Ary^/^, respectively. It 
is the first angle that tends to zero faster. If we take this an¬ 
gle exactly zero, as we have done in this paper, the scaling 
behaviour of the maximum violation will not change, but 
its value will be smaller by a factor approaching 6.25 near 
rjcrit, and by a factor of 1.33 at ry = 1 (see Eig. ([T]i). If we 
take the basis used at 12^ . given by the |000) and the \ W) 
states, the threshold efficiency remains 0.6, but near rjcrit 
we get much smaller violations; it will scale as the fourth 
power of A?y. This time 4>2 goes to zero faster than $ 1 . If 
we take $2 = 0, it will hardly affect the violation near rjcrit. 
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settings 

T^crit 

equation 

222 

0.6 


223 

0.6 

Gl 

333 

0.51678 


444 

0.5 

(GUI 


TABLE II. Table for critical detection efficiencies using symmet¬ 
ric 3-qubit states. The numbers in brackets refer to the Bell in¬ 
equalities. In the first column ijk refers to number of i,j, and k 
settings for the parties Alice, Bob, and Cecil. Detection efficiency 
thresholds are rounded up to five digits. 


while it will reduce it by about 30% at 77 = 1. The result is 
shown in Fig. ([T]) We may conclude that for this inequality 
our solution is much closer to the optimal arrangement than 
the one of Larsson et al. ll^ . However, their approach may 
directly be generalized to a larger number of parties. 


B. Symmetric state, m = 3 

The independent Bell coefficients we got for to = 3 are: 


= -2 

= -1 

Mill 

= 4 

iV(Zii2 — 

— 1 
ivjiia — 1 

= -2 

Mill 

= 1 

II 

CO CO 

— 1 
■^^^223 ~ ^ 

^2^33 = A 





For this Bell inequality Vent = (19 + v^)/96 
0.516776 (see Appendix), significantly smaller than the 0.6 
value we got with the \ W) state for to = 3. 


C. Symmetric state, m > 4 


For TO = 4 the coefficients are: 


Mlf 

= -2 

Mlf = -2 

-^'^-^112 

= 2 

Mill = 

-2 

M1^2 

= 2 

Mill = 1 

-^^^124 

= 1 

Mill = 

-2 


= 1 

— —2 

-^^^234 

= 1 

mI% = 

-2 

-^^'^334 

= 2 

-^'^344 ~ 




(28) 


In the Appendix we show that Vent is exactly 1/2 for this 
inequality. We have tried to = 5 and to = 6, but we have 
got no improvement, so for three participants we could not 
find a Bell inequality for which the critical efficiency goes 
below 1/2. 

We summarized critical detection efficiencies we found 
in this paper for the symmetric 3-qubit states in Table HU 
V. SUMMARY 

We have shown that the required detection efficiencies to 
demonstrate a loophole-free Bell violation can be signifi¬ 
cantly lowered if three parties are involved (instead of the 
usual two-party scenario). Before, no practical three-party 
Bell tests featuring efficiencies lower than 60% were known 
to the best of our knowledge. This value has been attained 
by Larsson and Semitocolos in 2001 in a three-party two- 
setting Bell scenario ll^ . We beat this limit using a W state 
and three measurements per party. Moreover, for 8 settings 
we reach the value of 50.13%. On the other hand, using a 
coherent mixture of the W state with a product state |111) 
allows us to obtain Vent = 50% even with 4 settings. We 
conjecture that Vent = 50% cannot be beaten in either way. 

It is left as an open question if one of our inequalities 
could be generalized beyond three parties similarly to the 
family of Bell inequalities by Larsson et al. 1^ . 
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Appendix A: Detailed calculation of critical detection 
efficiencies 

In the Appendix we calculate the critical detector effi¬ 
ciencies for the Bell inequalities given in the main text. Ac¬ 
cording to Eqs. mm, and taking into account the permu- 
tational symmetry of the states considered, the matrix ele¬ 
ment of the effective Bell operator may be written as: 

m m 

{-tplMnlit) ® (g) /|V') + 

j=l i=j 
m m m 

® Bj (g) C'fcIV'), 

k—1j=k i—j 

(Al) 

where is the number of permutations of indices i, j, 
and k, that is itijk = 6, if all three are different, = 3 if 
two indices agree, and Ttijk = 1 if i = j = fe. 

The condition for the violation of the Bell inequality by 
the results of the measurements performed on the \ W) state 
is: 

{W\M^\W) > 0, (A2) 

and the values of the matrix elements necessary to evalu¬ 
ate {W\J^ri\W) for small measurement angles are given by 
Eqs. (116117b . It makes the calculations simpler if we notice 
that these matrix elements do not change if we reverse the 
signs of the measurement angles concerned simultaneously. 
Also, if one of the measurement angles is $ and another one 
is — then the three particles matrix element will not de¬ 
pend on the third angle. These statements are also true for 
the matrix elements shown in Eq. (l22T i in the limit of small 
angles, which we will need when we consider the state de¬ 
fined by Eq. dZTli . 
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We have already shown that for the m = 3 inequality 
given by Eq. (O rjcrit = 0.6. 

Now let us consider the m = 4 case given by Eq. (HSJ. 
The measurement angles to be taken now are = (j)iX, 
with (/)i = — (/>2 = 1 and = —<^4 = A. By using 
Eqs. (IA2IA1I16I17I) . straightforward calculation leads us to: 

- 3 - A^ + 77[6 - 3(1 - 2Xf] > 0 (A3) 

for the condition of the quantum violation. Here we have 
simplified the expression by a factor of j:^? 7^/(48 • 12). At 
'H = Vcrit the l.h.s. of the equation is zero, therefore rjcrit = 
3(1 + 4A — 4A^)/(3 + A"^). It has its minimum value if A 
satisfies 2A® — 3A^ — A^ — 6 A + 3 = 0. The appropriate 
root calculated numerically is A = 0.466715, which leads 
to r]crit = 0.509036. 

Eor the the m = 6 case shown in Eq. (fTOl l we can follow 
the same steps as above. Now the measurement angles are 
given by 4>i = -4>2 = l,4>3 = -4>i = M and 4’b = -<('6 = 
V . With these angles we get for the condition of quantum 
violation, after a simplification by a factorof 77^/(48-36): 

—3(1 + + v'^) + 77[6 + 6/r"^ + 4i/^ — 3(1 — 2p)^— 

3(p2 - 2^jivf - 3(p - 7^ - pi/)2] > 0. (A4) 

Again, at 77 = rjcrit the l.h.s. of the equation is zero, and 
we must choose the parameters p and v such that rjcrit is 
minimal. We get three equations for the three unknown 
values, and if we solve those equations numerically we get 
ji = 0.495815, V = 0.295435, and rjcrit = 0.502417. 

Eor inequality with m = 8 given by Eq. (l20l i the expres¬ 
sion corresponding to Eq. dM is: 

—3(1 + r^) + 77(6 + 6 p^ + 6 cr"^ + 4 t^— 

3(1 - 2af - 3(p2 - 2paf - 3{a^ - 2 ctt)^- 
3{p — a — pa)^ — 3{pa — pr — ar)^] > 0. (A5) 


Here we have followed the same steps as for m = 6 taking 
measurement angles (jji = —(j )2 = 1, (/a = —(j)^ = p, 
05 = —06 = a and 07 = —07 = t. Erom the equation 
we get numerically rjcrit = 0.501338 with p = 0.498442, 
cr = 0.306395 and r = 0.169989. 

Now let the state be the one shown in Eq (|2TI) . Erom 
Eq (|26]) . if we choose the optimal mixing angle, the condi¬ 
tion for quantum violation is: 


{W\Mr>\W) 


(1E|7W^|111) 

(111|7W0111) 


(A6) 


The matrix elements of the Bell operator may be calculated 
from Eq (lAll i. which is also valid if the state vectors are 
different in the bra and the ket positions, provided both are 


permutationally symmetric. The matrix elements of the two 
and three particle operators appearin g in the r.h.s. of the 
equation are given in Eqs. (I16I17I22I) . We are concerned 
with the small angles limit. 

Eirst, let us take the m = 2 inequality of Eq (fT3]) . 
With the choice of 0i = 0 and 02 = 1, we 

get {W\Mr,\W) = -3773x^/48, (1E|7W^|111) = 

and (lll|Ai^|lll) = —irj^ + 2rj^ for the 
matrix elements of the effective Bell operator. By substitut¬ 
ing these values into Eq. (IA6I) . and taking into account that 
the l.h.s. of the equation is zero at 77 = rjcrit, h is easy to 
see that rjcrit = 3/5 = 0.6. 

We may take the same steps for m = 3. The inequality is 
shown by Eq. (|27]i, and the measurement angles are given 
by 01 = 0, 02 = 1 and 03 = —1. Then the matrix el¬ 
ements of the effective Bell operator are (W\J^r)\W) = 
—677^x^/48, (1E|A4^|111) = —y/3rj'^x^{2rj — 1/2) and 
(lll|A 4 ^|lll) = — 77^(12 — 577 ). Then the condition that 
the l.h.s. of Eq. (IA6b is zero at 77 = rjcrit leads to equa¬ 
tion A8rjlcit ~ ^^Vcrit — 3 = 0, whose appropriate root is 
Tjcrit = (19 -f v^)/96 Ri 0.516776. 

In the case of the 777 = 4 inequality of Eq. (l28l l the 
measurement angles are given by 0 i = 1, 02 = — 1, 
03 = A and 04 = —A. Erom these it follows that 
{W\Mr,\W) = rj^x^[-l -\^ + 77(1 -f 4A - 8A2 - + 

A^)]/4, {W\Mci\lll) = — •\/377^x^(1 -I- A3)(1 — 377 ) and 
(llljAf^|111) = — 2477^(1 — 77). If we substitute these 
values into Eq. (IA6I 1. we can get: 


8 ( 1 - 77 ) 



+p{v 



> 0 , 


(A7) 


where 


r = AX* + 8A3 -f 34A2 - 8A -f 7 
= 5A^ -f 2(A -f 1)'‘ -f 6A2 -f 4(2A - 1)^ -f 1 > 0 


p = 5A'‘ -f 6A2 -f 5 > 0 
(A2 -f 4A- 1)2 

0^- —4— 


(AS) 


In Eq. (IA7b the prefactor is positive for 0 < 77 < 1. As 
r is strictly positive, the inequality is satisfied above the 
upper root of the second order expression. Below that the 
expression is negative for all 77 > 0, as one can easily see. 
Therefore, we get the critical efficiency as {rjcrit — 1/2) = 
{y/p^ + Arq — p)/2r. As 7- > 0, p > 0 and q > 0, the 
smallest possible value the r.h.s. may take is zero, when we 
choose A such that <7 = 0, that is A = — 2 ± y/5. With this 
optimal choice rjcrit = 1/2. 












